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Abstract. We examine the convective core helium burning phase of intermediate mass stars, and investigate the role of coupling 
of nuclear burning and mixing on the extension and duration of the blue loop phase. We compare the theoretical scenario with 
the distribution of stars in the colour-magnitude (CM) diagram of the Large Magellanic Cloud (LMC) cluster NGC1866, whose 
largely populated clump of He-burning stars is equally populated in the blue and red side. We compare the distributions expected 
by adopting either a diffusive scheme within the instability regions, in which nuclear burning and mixing are self-consistently 
coupled, or the traditional instantaneous mixing approximation. We analyze with particular care the sensitivity of the results 
to: a) the e-folding distance with which the velocity of convective eddies decays outside the formal border of the convective 
regions fixed by the Schwarzschild criterion; b) the convective model adopted to evaluate the temperature gradient; c) the rate 
of reaction l2 C + a — » 16 O. Models not including convective overshoot are also commented. 
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1. Introduction 

The young stellar cluster NGC1866 in the Large Magellanic 
Cloud has been extensively studied in the last decades because 
it provides a unique opportunity to study the physics of the in- 
teriors of intermediate mass stars, particularly for what con- 
cerns the extension of the convective cores during the phases 
of both hydrogen and helium burning. It is young enough 
(Age <~ 200Myr) to have turn-off masses of approximately 
M ~ 4.5M Q , which are well known to develop extensive con- 
vective cores during their main sequence phase. At the same 
time the data available show an extraordinarily well populated 
clump region (Testa et al. 1999; Brocato et al. 2003), with 
~ 100 He-burning stars, and allow to compare the relative dura- 
tion of the permanence of He-burning stars in the red and blue 
parts of the CM diagram. 

The observational data have been used up to now primar- 
ily to find out informations regarding the convective central re- 
gions during the phases of central burning (Becker & Mathews 
1983; Brocato & Castellani 1988; Chiosi et al. 1989; Brocato et 
al. 1994; Barmina et al. 2002; Brocato et al. 2003), and particu- 
larly on the extension of the so-called overshooting region, i.e. 
the zone behind the formal border fixed by the Schwarzschild 
criterion in which the convective eddies are allowed to travel 
due to inertia despite the opposition of the buoyancy forces 
(Cloutman & Whitaker 1980; Xiong 1980; Canuto et al. 1996; 
Bressan etal. 1981). 
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The computations of detailed luminosity functions have 
provided a powerful tool to infer the relative duration of the 
phases of helium and hydrogen burning, which is directly 
correlated to the extension of the convective core during the 
MS phase, and the amount of extra-mixing which must be 
taken into account. Unfortunately, at the moment no common 
agreement has yet been reached, and the photometric study of 
NGC1866 has been presented by different groups either as a 
clue of a considerable amount of extra-mixing beyond the bor- 
der of the instability regions (Barmina et al. 2002), or as a clear 
evidence that such a region is so narrow that overshooting can 
be completely neglected in the evolutionary calculations (Testa 
et al. 1999; Brocato et al. 2003). To understand the level of un- 
certainty, suffice it to say that the differences in the ages of the 
cluster quoted in the literature amount to ~ 50%, ranging from 
~ lOOMyrupto ~ 150Myr. 

The idea behind this paper is not to provide another esti- 
mate of the age of the cluster through a comparison either be- 
tween the theoretical isochrones and the observed CM diagram 
or between the luminosity functions; rather, we focus our at- 
tention on the observed distribution of He burning stars in the 
CM diagram, in an attempt to reproduce the observed ratio be- 
tween blue and red clump stars (B/R ~ 1). We also discuss and 
compare the results obtained by varying the micro and macro- 
physics inputs adopted to calculate the evolutionary sequences. 

Our main goal is to check the sensitivity of the global du- 
ration of the helium burning phase and of the time spent in 
the blue region of the CM diagram on the way nuclear burning 
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and mixing of chemicals are coupled in the central convective 
regions during the two major phases of nuclear burning. We 
compare the results obtained by assuming a diffusive approach, 
in which the two processes are self-consistently coupled, with 
the instantaneous mixing approximation, still largely adopted 
in most of evolutionary codes. 

We use the observed distribution of clump stars in 
NGC1866 as a test of the diffusive scheme for the masses in- 
volved; we postpone to a forthcoming paper a more general 
discussion concerning a more extended mass range. 

We complete our study by testing the sensitivity of the re- 
sults also on the convective model adopted, our code allowing 
to calculate the temperature gradient within convectively unsta- 
ble regions through the FST or the MLT models (see Sect. 2.1), 
the rate of the reaction l2 C + a — > 16 O and the e-folding dis- 
tance of the convective velocity from the formal border of the 
convective core, fixed by the Schwarzschild criterion. 

A general description of the code ATON used to compute 
the evolutionary sequences is given in Sect.2. In Sect.3 we de- 
scribe the main evolutionary phases of the typical intermediate 
mass star which is likely to populate the clump of NGC1866. 
A comparison between the results obtained by varying the 
scheme for mixing is presented in Sect.4. Sects. 5 to 7 are fo- 
cused on the effects of changing, respectively, the e-folding 
distance of convective velocities decay within radiatively sta- 
ble regions, the convective model, and the cross-section of the 
reaction 12 C + a -> 16 O. 

2. The evolutionary models 

The stellar models calculated in this paper have been obtained 
by using the ATON2.0 code described by Ventura et al. (1998). 
The interested reader may find on the afore mentioned paper 
a detailed description of the numerical structure of the code, 
and of the macro and micro-physics which is used to simulate 
the stellar evolutions. Here we briefly recall the most important 
inputs adopted. 

2.1. The treatment of convection 

The calculation of the temperature gradient in stellar regions 
unstable to convection has been traditionally performed by 
adopting the Mixing Length Theory (MLT) (Vitense 1953; 
Bohm-Vitense 1958), which is based on the simplified assump- 
tion that the convective flux is carried by eddies which first 
form and travel for an assumed distance (the mixing length), 
then lose their identity, transferring to the surrounding layers 
their energy content. Both the dimension of the afore men- 
tioned eddies and the mixing length are assumed to be directly 
proportional to the distance of exponential decay of pressure 
(H p ), with a constant of proportionality omlt which is com- 
monly calculated in order to reproduce the evolution of the Sun, 
i.e. to fit the solar radius and luminosity at the present age of 
the Sun: we find omlt = 1-7. 

Alternatively, still in the framework of local convective 
models, in the last decade a new prescription was presented: 
the Full Spectrum of Turbulence (FST) model (Canuto & 
Mazzitelli 1991; Canuto et al. 1996). The main differences with 



respect to MLT are that all the spectrum of eddies' dimensions 
is taken into account, and the mixing length is assumed to be 
simply the distance to the closest convective border (/ = z). 
This latter model provides convective fluxes which are approx- 
imately ~ 10 times larger than the MLT ones in regions of high 
convective efficiency, and sensibly lower in stellar zones where 
the efficiency of convection is very poor, i.e. in proximity of the 
formal border or in the sub atmospheric layers. 

The code ATON2.0 can be easily managed in order to use 
either of the convective models. 

2.2. Nuclear network 

With respect to Ventura et al. (1998) the nuclear 
network has been extended and now includes 30 
chemicals:i/,D, 3 #e, 4 #e, 7 fie, 1 Li, 12 C, n C, H N, l5 N, lb O, 17 0, 
^O, m F, 19 F, 20 Ne, 21 Ne, 22 Ne 22 Na 23 Na 24 Mg 25 Mg 26 Mg, 
26 AZ, 27 A/, 28 S/, 29 S/, 30 S/, 3I P,«. 

All the main reactions of the proton-proton chain, CNO, 
Ne-Na and Mg-Al cycles, and of He burning have been consid- 
ered, for a total of 64 reactions. The relevant cross-sections are 
taken either by Caughlan & Fowler (1988) or by Angulo et al 
(1999). 

2.3. Convective mixing and nuclear burning 

The extension and the borders of all the convective regions are 
established by the Schwartzschild criterion. 

The treatment of nuclear burning within convective regions 
is based traditionally on the instantaneous mixing approxima- 
tion: it is assumed that convection is so fast in mixing the chem- 
ical species within the whole instability region that this latter 
can be assumed to be always homogenized. From a numerical 
point of view this latter assumption provides a great simplifica- 
tion, because it is possible to calculate average chemical abun- 
dances and nuclear cross-sections, which are used for all the 
mesh points included in the convective region. 

Alternatively, it is necessary to solve for each element the 
diffusion equation (Cloutman & Eoll 1976): 

(f)=(fL^[« 2D Sl a) 

stating mass conservation of element i. The diffusion coeffi- 
cient D is 

D = \6n 2 r A p 2 T- 1 (2) 

The turbulent diffusion time scale, t, is related to the correla- 
tion between density and velocity according to: 

<Pi<>=-r d $ (3) 

The knowledge of the second order momentum in eq.(0 re- 
quires the solution of the Navier-Stokes equations, which is not 
available. Thus, it is customary to rely on a local approximation 
for D, that is 

D = \ul (4) 
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where u is the convective velocity and / is the convective scale 
length. 

Careful details of the way of solving eq. Q can be found in 
Ventura et al. (1998) (see, in particular, the Appendix and Sect. 
2.3). 

A diffusive approach was used by Deng et al.(1996a;b) and 
Salasnich et al.(1999) to study the evolution of massive stars of 
solar and LMC metallicity during their phases of hydrogen and 
helium burning. 

The two schemes to treat nuclear burning within convec- 
tive regions need also two different approaches to deal with 
overshooting. In the instantaneous mixing approximation it is 
simply assumed, for the purpose of mixing, that the convective 
eddies travel up to a distance which is some fractions of H p 
away from the formal convective border. The numerical treat- 
ment is unchanged, the only difference being the larger exten- 
sion of the fully homogenized region. 

We thus have lov = <xH p . According to Maeder & Meynet 
(1991), a value of a — 0.25 is needed to fit the main se- 
quences of young clusters; for older clusters the observed data 
are hardly fit with values of a exceeding 0.2 (Stothers & Chin 
1992). 

Within the diffusive framework it is necessary to spec- 
ify the way with which convective velocities decay outside 
the convective boundaries (Deng at al. 1996a,b; Herwig et al. 
1997; Ventura et al. 1998). In agreement with Xiong (1985), 
Grossman (1996) and on the basis of numerical simulations by 
Freytag et al. (1996), we assume that convective velocities de- 
cay exponentially out of the formal convective boundary as: 

u = u » ex p ± (ih ln K) & 

where Ub and Pf, are, respectively, turbulent velocity and pres- 
sure at the convective boundary, P is the local pressure, ( a free 
parameter connected with the e-folding distance of the decay, 
and f thick is the thickness of the convective regions in fractions 
of Hp. The chemical abundances are still calculated by eq.Q, 
which is applied in a slightly more extended region with respect 
to the borders fixed by the Schwarzschild criterion. 

a and £ are the two parameters connected with the extra- 
mixing in both schemes, though their physical meaning is com- 
pletely different. In terms of duration of the main sequence 
phase for intermediate mass stars, we may say that f = 0.02 
leads to results similar to a — 0.18, and ( = 0.03 to a ~ 0.25. 

2.4. Chemical composition 

We calculated several sets of stellar tracks by varying the 
scheme adopted for chemical mixing (instantaneous or diffu- 
sive), the convective model (FST and MLT), the parameter for 
the exponential decay of velocities = 0.02 and £ = 0.03) 
and the rate of the reaction 12 C + a — > 16 O. Both for the 
instantaneous and the diffusive case a set of evolutionary se- 
quences where the extension of the mixed region is tightly fixed 
according to the Schwarzschild criterion was calculated. The 
adopted chemical composition for all the tracks discussed here 
is Z = 0.008 and Y = 0.25 (Brocato et al. 2003; Barmina et al. 
2002). 




B-V B-V 

Fig. 1. The theoretical isochrones superimposed on the obser- 
vational data for LMC cluster NGC1866. Left: isochrones cal- 
culated with diffusive mixing and £ = 0.02. Right: instanta- 
neous mixing, a = 0.18. 



3. Some evolutionary remarks 

For all the sets of tracks considered we calculated theoretical 
isochrones for ages in the range 100 - 200Myr following the 
scheme suggested by Pols et al.(1998). The conversions from 
the theoretical to the observational plane have been performed 
by means of the transformations by Castelli et al.(1997). We 
first tried to select for each set of tracks the most appropriate 
parameters (age, distance modulus, reddening) allowing a de- 
cent fit of the observed CM diagram. First, we checked the pos- 
sibility of discriminating between the two schemes for mixing, 
looking for differences in the quality of fitting obtained with 
both prescriptions. 

Fig-Eshcws the comparison between the observed CM di- 
agram (Testa et al. 1999) with the theoretical isochrones ob- 
tained by assuming a diffusive scheme for mixing with ( = 0.02 
(diffl)2 models), or an instantaneous prescription with an over- 
shooting distance of lov = 0.18H P (istl8 models). The con- 
vective model adopted is FST, and the nuclear network is from 
Caughlan & Fowler (1988). 

In both cases we see that a decent fit of both the extension 
of the main sequence and the luminosity and extension of the 
clump is obtained with an age of 150Myr, a distance modulus 
of (m - M) = 18.6 and a reddening E(B -V) = 0.05. 

This preliminary approach seems to rule out the possibility 
of selecting different mixing modalities simply on the basis of 
CM fitting. A deeper analysis is needed. Before entering into 
the details of this comparison, we briefly examine the physical 
properties of the stars which are supposed to populate both the 
turn-off and the clump region within the context of our sim- 
ulations. This will help us to better understand the results we 
obtain. 
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M = 4.5M S 
Z=0.008; Y=0.25 
<-=0.02 





Fig. 2. The evolution on the theoretical HR diagram of a model 
of M = 4.5M calculated with diffusive mixing = 0.02). 



The typical mass belonging to the clump of fig.^is 4.5M . 
The track on the HR diagram is shown in fig. |2j the bottom 
panel shows in more details the clump region. In the early 
phases of hydrogen-burning the star develops a central convec- 
tive region of ~ 1.1M with a radial extension of ~ 0AR Q . As 
CNO burning proceeds, the central regions contract and heat, 
the convective core shrinks in mass (it drops to ~ O.5M when 
the abundance of hydrogen is 0.25 in the centre) while its radius 
decreases slowly. Consequently, the pressure scale height at the 
border of the central convective region diminishes, and the ex- 
tension of the extra-mixing region, directly connected to H p , is 
also reduced. The whole phase of H-burning lasts ~ 130Myr, 
after which the star expands. Due to the lower and lower effec- 
tive temperatures convection penetrates deeply towards internal 
regions which were previously touched by nuclear burning, in 
what is commonly known as the first dredge-up. The base of the 
external envelope reaches a layer which is ~ 1M in mass away 
from the centre, while the stellar luminosity, due essentially to 
CNO burning in a shell, rises up to log{L/L Q ) ~ 3.2, until 3a 
reactions start in the centre of the star (point A in the bottom 
panel of fig. 01. At the ignition of helium in the centre the en- 
ergy release is in excess with respect to the rate of energy which 
can be carried outwards. This triggers a rapid expansion of the 
central core, which has the effect of diminishing the strength of 
the CNO shell burning, so that the luminosity decreases. The 
stellar layers from the shell to the surface contract, convection 
recedes, the track moves downwards in the HR diagram until 
the CNO shell, after ~ 10 7 yrs since the beginning of helium 
burning, meets the chemical discontinuity left behind by the 
first dredge-up (point B). This determines a rapid increase of 
the CNO (and global) luminosity, so that the star continues to 
burn helium in the red part of the HR diagram until surface con- 
vection disappears (point C). Then a readjustment of the struc- 
ture causes a rapid contraction, and the tracks moves to the 



Fig. 3. The evolution with time of the internal structure of the 
same model of fig.|2]in terms of the coefficient of generation 
of nuclear energy (top) and the hydrogen profile (bottom) The 
numbers in the upper panel indicate value of the corresponding 
stellar luminosity, while characters B and D corresponds to the 
B and D evolutionary phases evidentiated in fig. [3 We note in 
the bottom panel the chemical discontinuity at M/M Q ~ 0.95. 



blue part of the diagram until the CNO shell evolves through 
the discontinuity (point D). The star remains in the blue part 
of the clump region until the helium abundance in the central 
regions drops to Y ~ 0.2, after which the rate of nuclear energy 
release drops and the core starts contracting; the layers from 
the CNO shell to the surface expand, and the tracks moves to 
the red. In the diff02 model the track stays in the bluest region 
of the clump for ~ 7Myr. The whole phase of helium burning 
lasts approximately ~ 20Myr. 

Fig. shows the variation with time of the coefficient of 
generation of nuclear energy (top panel) and the hydrogen pro- 
file (bottom) at some stages intermediate between points B and 
D of fig. |2 We note a rapid decrease of the energy release un- 
til the outer part of the shell reaches the chemical discontinu- 
ity (evidentiated by the sudden increase of the hydrogen abun- 
dance at M/M Q ~ 0.95) (point B); point D corresponds to the 
phase when the peak of the shell reaches the quoted disconti- 
nuity. 

We stress here the different behavior of the convective cen- 
tral regions of the star with respect to the preceding phases of 
hydrogen burning. During H-burning, the extension in mass of 
the core tends to increase as nuclear burning proceeds. The sec- 
ond important difference is that the ignition of 3a reactions pro- 
vokes a rapid expansion of the central regions, so that while the 
temperature increases, pressure and density decrease: this has 
the effect of increasing the value of H p at the border of this re- 
gion, with the consequence that more material external to the 
formal convective region is carried inwards. 
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This picture outlines which are the main factors driving the 
evolution of the star during the He-burning phase: 

- The excursion of the track to the blue is connected to the 
disappearance of the convective envelope. 

- The bluest point of the track is reached when the whole 
CNO shell has evolved through the chemical discontinuity 
left behind by the first dredge-up: we may thus expect that 
the extension and the morphology of the blue loop depend 
on the mass coordinate of the innermost point reached by 
external convection during the dredge-up and on the chem- 
ical stratification left. 

- The return of the track to the red is connected with the con- 
traction of the internal region just before helium exhaus- 
tion. The duration of the permanence of the star in the blue 
is therefore directly dependent on the amount of helium still 
present in the core when the star reaches the bluest point 
of the loop. The way with which central regions are refur- 
bished of helium are expected to play a fundamental role in 
determining the time spent in the blue. 

4. Instantaneous vs. diffusive mixing 

Fig.[2shows that a decent fit of the HR diagram can be easily 
obtained both with an instantaneous and a diffusive approach 
to deal with burning of chemicals in the regions unstable to 
convection, because no clear difference can be seen in the the- 
oretical isochrones built with the diff02 and istl8 models. 

To be able to distinguish between the results of the two 
schemes we have to build a synthetic CM diagram, and to com- 
pare it with the observed CM diagram by Testa et al. (1999). 
Our main goal is to reproduce the distribution of stars in the 
clump, i.e. to obtain the observed ratio of B/R ~ 1 between the 
blue and the red stars of this region. 

To build the synthetics, we developed a simple algorithm 
which populates a given isochrone with a random distribu- 
tion of masses, chosen to be available for the isochrone of 
that given age. For each extracted mass, values of LogL and 
LogTe, as well as magnitudes and colours in the observational 
planes, were obtained by means of linear interpolation from the 
points of the same isochrone. This simple way of populating an 
isochrone gives us important indications regarding the relative 
populations of the various evolutive phases, allowing us to de- 
rive the expected B/R ratio for each of our tests. 

This ratio was calculated by counting the number of stars 
within appropriate boxes in the CM diagram chosen on the 
basis of the observational distribution of He-burning stars in 
NGC1866. Specifically, we retained as blue stars those whose 
(B - V) and V fall within the range, respectively, of 0.5 < 
(B - V) < 0.9 and 15.6 < V < 16.3; we considered as red 
stars, those with 1 < (B —V) < 1.5 and 15.5 < V < 17.1. We 
deliberately chose not overlapping boxes in order to prevent 
small fluctuations in colour to substantially affect our statistics. 

Statistical fluctuations related to the random distribution of 
star masses are taken in account in the following way: for each 
choice of parameters, we performed ten simulations varying 
only the seeds for random number generation: after that, the 




0.5 1 1.5 0.5 1 1.5 

B-V B-V 



Fig. 4. Synthetic CMD computed by adopting diffusive stellar 
models with ( = 0.02 (left) and instantaneous models with a = 
0.18 (right). In the lower part we show the hystohistogram of 
the distribution in colour of the He-burning stars; we note the 
presence of a blue population with (B - V) in the range 0.4-0.5 
in the diff02 case. 



computation of B/R ratio is made on the averaged numbers of 
blue and red, according to the definition given above. 

Note that we did not include the age spread in the follow- 
ing simulations. However, preliminary computations show that 
a reasonable age spread (~ 10% of the age of the cluster) can- 
not alter significantly the derived B/R ratio (in agreement with 
Barmina et al. 2002, see their Fig. 15). 

Fig.|4]shows the numerical simulations obtained for the two 
sets of models: we can see that in the diff02 case the blue part is 
much more populated, with a ratio B/R ~ 60%, to be compared 
to the B/R ~ 25% value corresponding to the istl8 models. 
This can be clearly seen in the corresponding histogram shown 
in the lower part of fig.E] Such a steep dependence of the time 
spent in the blue and red side of the CM diagram of He-burning 
intermediate mass stars on the treatment of mixing within the 
convective core was already discussed in D' Antona (2002). 

To understand the reason of this difference we compared 
carefully the evolutions of a 4.5M model of the diff02 set with 
a istl8 model of 4.4M , which are the typical masses staying 
in the clump at the age of 150Myr for the two sets of model 
discussed. The two tracks are very similar during the MS evo- 
lution (see fig. [5}, and have approximately the same luminosity 
during the helium burning phase. The duration of the H-burning 
phase is the same (~ 130Myr) for both models. 

A detailed comparison between the two evolutions can be 
seen in the three panels of fig. [6] In the abscissa we report 
the time starting from the beginning of hydrogen exhaustion 
in the centre; to have an idea of the extension of the loop and 
the duration of the blue phase, we show in the lower panel the 
variation of the effective temperature. We find in both cases a 
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Fig. 5. The comparison of the tracks follower in the HR dia- 
gram by two models calculated with two different prescription 
for mixing. The selected masses populate the clump region ac- 
cording to our simulations. 
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Fig. 6. The variation with time of some quantities related to the 
evolution of the same models shown in fig|5] In the abscissa we 
report the time counted from the hydrogen exhaustion in the 
stellar centre. Upper: Central helium; Middle: Central density. 
Lower: Effective temperature. 

rapid increase of T e fj with time in conjunction with the disap- 
pearance of surface convection. The diffusive track stays in the 
blue longer (tbi ue ~ 7Myr), while the istl8 track, as soon as it 
reaches the bluest point, moves again to the red (tbi ue ~ 3Myr). 

The reason of such a different behavior can be detected in 
the first and second panel of the same figure. In the upper panel 
we see that in the istl8 model helium is burnt faster, so that 



when the bluest point is reached the central helium abundance 
has dropped to ~ 0.2; with these low helium abundances the 
two nuclear reactions involved (3a and 12 C + a — > 16 O) can- 
not supply the whole energy demand, so that the contraction 
of the core, triggering the expansion of the external layers and 
the return of the track to the red, starts soon. In the diffusive 
model the helium burning is slower, so that at the bluest point 
there is still 40% of helium in the core, thus delaying the be- 
ginning of core contraction. This allows the track to stay in the 
blue region for longer, and the relative duration of this phase is 
approximately doubled. 

The results are in full agreement with the detailed com- 
parison of the H and He-burning times obtained with the two 
schemes for mixing for the most massive intermediate mass 
stars presented in Salasnich et al.(1999) (see their fig. 5). 

These results indicate that the two schemes for mixing can 
lead to deeply different results, particularly in terms of the rel- 
ative duration of the various evolutionary phases. The diffusive 
scheme allows an exponential decay of velocities beyond the 
border of convection, with the consequence that the zone which 
is mixed is larger: yet, we must consider that the time-scales for 
mixing in the most external layers are extremely large, so that 
the refurbishment of fresh helium into the hottest central lay- 
ers is slow. On the contrary, in the instantaneous case a smaller 
region is involved in mixing, but the efficiency of convection 
in this latter case is assumed to be so high to homogenize it 
completely: the process of carrying helium to the centre is thus 
more efficient, the rates of both 3a and 12 C + a — > 16 O reac- 
tions are increased, and the whole process of helium burning is 
faster. 

It is not surprising that the differences between the results 
obtained with the two schemes for mixing are more evident 
during the He-burning phase, since in this case the size of the 
convective core tends to increase; during the H-burning phase 
the shrinking of the convective core, precedently discussed, 
makes the size of the extra-mixing regions to be progressively 
decreased, so that the differences between the two treatments 
are smaller. 

The comparison between the results obtained with differ- 
ent schemes for mixing leads to the conclusion that interme- 
diate mass models computed with the diffusive approach, pro- 
vided the luminosity of turn-off and of the clump are the same, 
consume helium more slowly than the corresponding models 
computed with the instantaneous mixing approximation, with 
a consequent longer duration of the blue phase. 

Applied to the study of NGC1866, the usage of the diffu- 
sive approach allows to shorten the gap between theory and 
observations, rising the predicted B/R ratio to ~ 60% (to be 
compared to the expected value B/R ~ 25% on the basis of the 
istl8 models). 

4.1. The no extra-mixing case 

The same analysis to test the influence of the time-dependent 
algorithm for mixing was applied to two sets of tracks calcu- 
lated by assuming that no sort of extra-mixing takes place out 
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Fig. 7. Evolutive tracks of istOO models with masses in the 
range 3.5 - 5M in the theoretical plane. We note the very nar- 
row extension of the loop for the most massive models shown. 



of the border fixed by the Schwarzschild criterion (diffOO and 
istOO models). 

Fig. shows the tracks followed in the HR diagram by 
models with masses in the range 3.5 - 5M , calculated with the 
instantaneous mixing approximation. We note that the bluest 
point of the track gets hotter as the stellar mass increases, but 
for models more massive than ~ 4.3M the track substantially 
stays in the red during the whole He-burning phase; this behav- 
ior is found for all models less massive than ~ 6.5M . These 
results are in full agreement with the most recent findings con- 
cerning the evolution of intermediate mass stars of this metal- 
licity by Pietrinferni et al. (2004). 

This can be explained on the basis of the fact that in the 
more massive models helium consumption if faster (due to the 
higher temperatures) and convection is always well extended: 
when the abundance of helium drops to Y ~ 0.2, and the core 
contraction begins, the surface convection is still alive, so that 
the star never begins the rapid excursion to the blue which is 
typical of models with radiative envelopes. Therefore, we find 
no blue clump stars for ages t < 140Myr. 

Disregarding the presence or not of the blue loop, on the 
basis of the clump luminosity we obtain an age for NGC1866 
of ~ lOOMyr with a distance modulus of (m - M)o ~ 18.5, in 
agreement with Testa et al. (1999); for a lower distance mod- 
ulus ((;« - M)o ~ 18.25) we obtain a decent fit of the clump 
luminosity for an age of ~ 150Myr (in agreement with Brocato 
etal. 2003). 

The lack of stars burning helium in the blue seems to rule 
out ages < 140Myr, and in any case does not allow us to per- 
form the usual analysis based on the B/R ratio. For an age 
t = 150Myr we cannot fully reproduce the extension of the 
clump both in temperature and luminosity; in any case, the 
percentage of blue stars compared to the red would be low 



(B/R ~ 25%) compared to the observed value of ~ 1. This can 
be understood on the basis of the helium still present within 
the core when surface convection is extinguished. The lack of 
any extra-mixing out of the formal convective border triggers a 
rapid consumption of helium in the central region, so that when 
the excursion towards the blue begins the abundance of helium 
dropped to Y ~ .25 - .3: soon after the contraction of the central 
layers forces a quick return of the track to the red. 

It's interesting to note that the B/R ratio expected on the 
basis of the instantaneous models calculated with no extra- 
mixing is the same as in the istl8 models, precedently dis- 
cussed. Therefore, within the framework of instantaneous mod- 
els, it is not possible to derive informations concerning the ex- 
tension of the overshooting region on the basis of this kind of 
analysis only. 

If we compare the results obtained with the two schemes 
for mixing (diffOO and istOO models), we note that in these 
models calculated by assuming a complete absence of extra- 
mixing no meaningful differences have been found, both in the 
main sequence and in the helium burning phase. This can be 
explained on the basis of the fact that the time scale of convec- 
tion is very short with respect to those typical of the various 
reactions involved within the formally unstable region. For ex- 
ample, within the 4.3M istOO model, which would populate 
the clump of NGC1866 if no extra-mixing is considered, the 
time-scale of convection is always shorter than ~ 1 month for 
the whole duration of the He-burning phase. 

5. The role of f 

The discussions of the previous sections outlined the difference 
between the roles of the parameters £ and a in determining 
the amount of matter which is mixed and carried to the central 
regions during the major phases of central burning. 

( is connected with the e-folding distance of the velocity 
decay from the border of the convective zone, but even in case 
of very large values of £ material far from the border is hardly 
mixed to the centre, because the time-scale for mixing in that 
region would be comparable or even longer than the evolutive 
times. Hence, we expect that even large percentage variations 
of C, don't have dramatic consequences on the evolution of the 
star, unlike the effect of great changes of a. 

A preliminary study regarding the effects of varying the 
e-folding distance of the exponential decay of velocities out- 
side the border of the convective core is found in Salasnich et 
al.(1999), but it is mainly focused on the evolution of massive 
stars (M ~ 2OM ). A straight detailed comparison between our 
parameter £ and their quantity a\ is made difficult by the differ- 
ent way of simulating the exponential decay: we use pressure 
as independent variable, while they use the radius. 

We increased the value of £ up to 0.03, obtaining a nice 
fit of the CM diagram of NGC1866 with the same param- 
eters holding for £ = 0.02, but with a slightly larger age 
(160Myr). In this latter case the mass evolving in the clump 
is Mclump ~ 4.4M . The numerical simulations run for this 
set of tracks show (see fig.|8} that the percentage of blue stars 
is increased up to B/R ~ 75%, which is still lower than the 
observed value, but pointing in the right direction. 
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Fig. 8. Synthetic CMD computed by adopting diffusive stellar 
models with ( = 0.03 (right) compared with the observed CMD 
of NGC1866 by Testa et al. (1999). In the lower part of both 
panels we show the corresponding histogram related to the dis- 
tribution in colour of the He-burning stars. 



The increase of the time spent in the blue for larger values 
of £ can be understood again if we consider the fast consump- 
tion of helium within the central regions, and the quantity of 
helium still left when the track reaches the bluest point of the 
clump. 

To understand the differences between the ( = 0.02 and £ = 
0.03 sets we compare two models of masses M = 4.4M (case 
£ = 0.03) and M = 4.5M (£ = 0.02). Such a choice stems 



Fig. 9. Path followed by two models calculated with the dif- 
fusive scheme for mixing and different values of ( which are 
currently populating the clump region according to our simula- 
tions. 

from the small difference found between the masses populating 
the clump region in the best simulations obtained with the two 
sets of models. 

In fig.[9]we compare the tracks on the HR diagram; we note 
that both the luminosity of the turn-off and of the bluest point 
of the clump are almost the same. In the 6 panels of fig.llOlwe 
show in details the differences in terms of various chemical and 
physical quantities during the He-burning phase. 

In the left panels we report the evolution with time of the 
CNO burning shell in terms of the mass coordinate of the layer 
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corresponding to the maximum release of CNO nuclear output 
(left - upper panel) and of the maximum value of the coefficient 
for the generation of CNO nuclear energy (left - middle panel). 
The right panels shows the evolution of the central physical- 
chemical conditions, in terms of density, temperature, and he- 
lium abundance. 

In the f = 0.03 model the H-exhausted region is larger, 
and the CNO shell is more external (in mass); the difference 
is of the order of Am ~ O.O4M . The CNO burning shell 
in the £ = 0.03 model, before the ignition of 3a process, is 
consequently hotter and more efficient, so that the energy re- 
lease is larger, as we can see in the left - middle panel. The 
faster drop of CNO energy release which can be seen in the 
( - 0.03 model, is related to the conditions in the core. The 
right-middle panel of fig. [K)] evidentiates that the central re- 
gions of the f = 0.03 model are slightly hotter, so that helium 
is consumed faster (right-lower panel); this, in turn, triggers a 
more rapid expansion of the central regions, which can be de- 
tected in the steeper drop of central density in the right-upper 
panel, and a more rapid decrease of the nuclear energetic output 
due to CNO burning shell (see the left-middle panel). 

The evolution of the £ = 0.03 model, for what concerns 
the conditions of the external layers of the structure, is faster, 
so that the overall contraction leading the track to the blue in 
the HR diagram is anticipated by ~ 2Myr with respect to the 
( = 0.02 case. 

The same difference in the velocity of the evolution is not 
found in the central regions, because the larger temperatures 
within the interior of the ( = 0.03 model (which would lead to 
a more efficient nuclear burning and thus to a faster helium con- 
sumption) are partly balanced by the larger value of f, which 
allows to refurnish the core with more helium from the regions 
away from the formal border of the convective core to be trans- 
ported in the hot central regions. The earlier excursion to the 
blue of the £ = 0.03 track is not followed by a quicker return to 
the red (see the left - lower panel of fi g . 1 1 Oi so the percentage 
of time spent in the blue increases (~ 70%, to be compared to 
~ 50% for £ = 0.02). 

We also note that the bluest point of the ( = 0.03 model 
nicely fits the location on the HR diagram of the bluest Cepheid 
EV Set found in the open cluster NGC 6664 (Schmidt 1984). 

These results, together with the conclusions drawn at the 
end of the previous section, show that the diffusive model with 
£ = 0.03 leads to a better agreement between the theoretical 
framework and the observational scenario. Thus, we decided 
to adopt this set of models to perform the further analysis to 
test the influence of the convective model and of the relevant 
reaction rates. 

6. Can the CM diagram of NGC1866 help in 
discriminating between convective models? 

The determination of the temperature gradient (or, equivalently, 
of V = jj^y) within convective regions is one of the most 
disputed issues in modern stellar astrophysics (Ulrich 1976; 
Grossman 1996; Grossman & Narayan 1993; Xiong et al. 1997; 
Canuto 1997; Canuto & Dubovikov 1998). 



Age=160Myr; f=0.03 
FST 




Fig. 11. The comparison between two isochrones of the same 
age corresponding to models calculated with two different 
convective schemes. We note the different slope of the giant 
branches, the MLT isochrone becoming cooler and cooler for 
lower T e ff's. 

These difficulties stem from the absence of an adequate so- 
lution of the Navier - Stokes equations, so that the physical 
description of the instability regions still rely on local approx- 
imations, in which V is calculated on the basis of the values 
which some physical quantities assume in the same point where 
we want to know V. 

This approximation neglects one of the most important fea- 
tures of turbulence, i.e. non locality, but still it is the only op- 
portunity allowing the evaluation of V. 

Canuto & Mazzitelli (1991) constructed a local model for 
turbulent convection (FST model) which has been successfully 
employed in several astrophysical studies, like helioseismol- 
ogy (Basu & Antia 1994), red giant evolution (Stothers & 
Chin 1995; 1997), AGB evolution of intermediate mass stars 
(D'Antona & Mazzitelli 1996; Mazzitelli et al. 1999). Within 
this context we test the possibility of selecting between the 
two local convective models currently available, i.e. the Mixing 
Lenght Theory (MLT) and the FST model, on the basis of 
the comparison between the theoretical predictions and the ob- 
served distribution of NGC 1866 stars in the CM diagram. 

We therefore calculated a set of MLT models with a dif- 
fusive scheme for mixing; the parameter for the exponential 
decay of convective velocities was set to £ = 0.03, on the basis 
of the discussion presented in the previous section. 

In fig.[n we compare the isochrone giving the best fit with 
the £ = 0.03 FST models with the corresponding isochrone 
built by using the MLT models. The two lines are almost per- 
fectly superimposed during H-burning, while we note a differ- 
ent slope of the giant branch, the MLT line becoming progres- 
sively cooler for lower effective temperatures. The luminosity 
and the extension of the clump are very similar. 
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The numerical simulations (not shown here) show no mean- 
ingful differences in the distribution of stars in the clump, par- 
ticularly in terms of the B/R ratio; the only difference between 
the two sets of models is in the colour of the giant branch, 
which is redder in the MLT models. The sensitivity of the col- 
ors of giants to the convective model adopted has been explored 
by Stothers & Chin (1997) (see also the discussion in Asida 
2000), who found that a further tuning of the mixing length 
parameter with mass and effective temperature is needed to fit 
the observed slopes of red giant branches. We recall that in the 
present MLT computations the solar calibrated mixing length 
parameter was used for all the evolutions. 

On the theoretical point of view these results are indeed not 
surprising, if we consider the physical contexts where we ex- 
pect the convective model used to have an impact on the results 
obtained, in terms of temperature stratifications. Actually, in- 
side convective regions where the efficiency of convection is 
very large (i.e. almost all the energy flux is carried by con- 
vection), the gradient sticks to the adiabatic value, and this is 
almost independent of the description of convection adopted. 

In the case of NGC1866, if we consider the evolution of the 
masses currently in the clump according to our simulations (i.e. 
M ~ 4.5M ), we can select three different cases where the star 
develops wide convective regions: 

- The convective core during H-burning. 

- The convective central region which is formed along with 
the ignition of 3a processes. 

- The convective envelope which forms following the expan- 
sion of the external layers after hydrogen exhaustion. 

In the first case (core H-burning), the overadiabaticity (V - 
Vad) is on the average ~ 10~ 7 within the whole convective 
core in both models; the only difference can be noticed just in 
the very proximities of the border, where the FST (V - Vad) 
is seen to rise up to 10~ 5 , while in the MLT case we find 
(V - Vod) ~ 10~ 8 . This different behavior can be understood 
on the basis of the different mixing scales adopted in the two 
cases (Canuto & Mazzitelli 1991). In the FST model the choice 
A = z, the distance from the border, leads to a lower convec- 
tive efficiency, while the MLT, adopting a typical scale propor- 
tional to Hp, simulates a much larger efficiency of convection 
in these regions. Even taking into account these differences, 
we still have a very small degree of overadiabaticity in a region 
where the pressure drop is approximately 2 dex (log(P) ~ 17 in 
the centre, and log(P) ~ 15 at the border), so that in both cases 
the temperature gradient sticks to the adiabatic value. This can 
explain the similarity of the two isochrones in the MS region 
of the cluster. 

In the second case (core He-burning) we expect the differ- 
ences between the results provided by the two prescriptions for 
convection to be even less than during H-burning, because of 
the larger densities in the He-burning core, leading to a larger 
efficiency of convection. Indeed, we find in this case that along 
the whole convective core (V - V ac i) ~ 10~ 8 , with a maximum 
value of (V - V fl j) ~ 10~ 7 attained by the FST model in the 
regions next to the border. We may conclude that even dur- 
ing He-burning the temperature stratifications are very close, 



so that the differences in terms of duration of the whole He- 
burning phase are completely negligible. 

The third point quoted above (convective envelope) de- 
serves particular attention, because the outermost layers of the 
star, being cooler and less dense than the interiors, present val- 
ues of V which largely deviates from the adiabatic values. 

Within the convective envelope of the two models prece- 
dently discussed we note a mean value of overadiabaticity 
which is (V - V flrf ) ~ 3 x 10~ 5 in the FST case, and (V - V ad ) ~ 
3x1 0~ 4 in the MLT model: this difference is due to the larger ef- 
ficiency of convection of the FST model, which leads to larger 
convective fluxes for a given value of overadiabaticity (Canuto 
& Mazzitelli 1991). These values are considerable higher than 
those found within the central cores; this, together with the se- 
vere drop of pressure from the base of the external zone up to 
the surface (~ 5 — 6 dex) makes the differences in terms of 
temperature stratification more evident. 

In the very proximities of the surface the FST model has 
an overadiabaticity peak which is higher and narrower than the 
MLT one: this is again due to the MLT choice of the scale of 
mixing, which acts to simulate a convective efficiency which 
keeps sufficiently large even in regions where in reality a very 
low percentage of energy is carried by convection. Yet, these 
differences involve only the outermost layers, therefore do not 
have any incidence on the results obtained. 

We may conclude that in the MLT case the slope of the 
temperature drop from the base of convection to the surface is 
higher, so that the effective temperature is lower; this becomes 
more evident for deeper convective envelopes, so we may un- 
derstand the different slope of the giant branches which can be 
seen in fig. II II 

On making the comparison with the observed CM diagram 
of NGC1866, we don't find any appreciable difference in terms 
of the B/R ratio, for two main reasons: 

- During the evolution in the blue the only convective zone 
within the star is the He-burning core, for which we already 
stressed that the gradient is almost adiabatic: the evolution 
of the star in terms of evolutive time and colour is therefore 
independent of the convective model. 

- In the red region of the clump, we saw that the convective 
model directly influences the values of T e ff> but not of the 
thermodynamical quantities of the CNO burning shell; the 
time spent in the red is therefore not influenced by the con- 
vective model. 

This comparison between the observational scenario and 
the theoretical predictions for two convective models can only 
confirm the problems of fitting the colours of giants connected 
with the MLT model, but, since the influence upon the B/R ratio 
is very poor, it is not possible within this context to discriminate 
between convective models. 

7. The reaction 12 C + a -V 6 O. 

The evolution of a star in the early phases of He-burning is 
dominated by the 3a reactions, but as helium is consumed in 
the central regions the reaction 12 C + a — > 16 O becomes more 
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Fig. 12. The variation with temperature of the cross section (ex- 
pressed in cm^'mole -1 units) of the reaction n C + a — > 16 O 
according to the prescriptions given in Caughlan & Fowler 
(1988, CF88), Angulo et al. (1999, NACRE), and Kunz et al. 
(2002). 



and more important, eventually determining the duration of the 
last phases of helium burning and the residual chemical abun- 
dances left in the central regions of the star after helium ex- 
haustion. 

The cross-section cr\i c+a is uncertain by at least a factor 2, 
so we want to check the sensitivity of the expected distribution 
of He-burning stars in the CM diagram which may be due to 
changes in the adopted relationship <ri2 C+a (T). 

The influence of the quoted reaction on the times spent by 
intermediate mass stars in the helium burning phase, and par- 
ticularly on the blue region of the clump, has been explored 
by Imbriani et al.(2001), who compared the results obtained 
by adopted either the Caughlan & Fowler (1988, hereinafter 
CF88) or the Caughlan et al. (1985) rates, this latter release pro- 
viding approximately rates larger by a factor of 2. Their main 
findings, essentially, were that the tracks in the HR diagram 
are not dramatically altered by the choice of cr\i c+a , while the 
relative duration of the blue phase during helium burning is in- 
creased by ~ 10%. This result would thus increase the percent- 
age of blue stars in the clump region, helping to reconciliate 
the observational scenario for NGC1866 with theoretical pre- 
dictions. 

Here we want to compare the results obtained by adopt- 
ing the CF88, the NACRE (Angulo et al. 1999) and the Kunz 
et al. (2002) rates for the quoted reaction. We therefore calcu- 
lated the same stellar models giving the best fit of the observa- 
tional scenario (i.e. diffusive models with £ = 0.03) with the 
NACRE cr\i c+a (T) relationship, which provides values larger 
by a factor of ~ 2 over the whole range of temperatures of in- 
terest here (Angulo et al. 1999), and with the Kunz et al. (2002) 
prescription, which provides cross-sections sligtly lower than 
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Fig. 13. The path followed in the HR diagram by three models 
of mass M = 4AM G calculated with different cross-sections for 
the nuclear reaction l2 C + a — > 16 O. The solid track refer to the 
cross-section taken by Caughlan & Fowler (1988), the dotted 
line has been calculated with the NACRE rate by Angulo et al. 
(1999), and the dashed line with the rate given by Kunz et al. 
(2002). 
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Fig. 15. Variation of the central thermodynamical conditions of 
two of the three models presented in fig. during the early 
phases of helium burning. 

NACRE, and more sensitive to temperature changes in the 
range 10 8 < T(K) < 2 ■ 10 8 , which are those of interest here 
(see fig.ll2>. 

Fig. ^] shows the comparison between three tracks of a 
model of mass M = 4AM & , the typical mass populating the 
clump region (see Sect. 5), calculated with the three different 
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Fig. 14. The variation with time of various chemical and physical quantities related to the evolution of the same models shown in 

fig.na 



reaction rates. We focus only on the clump region, since no 
differences are found during the precedent phases. We see that 
the morphology and the extension of the loops are very similar, 
with the only difference that the NACRE model is slightly more 
luminous in the upper part, which corresponds to the phase of 
expansion of the star as the central helium abundance falls be- 
low Y ~ 0.1. 

To understand which results we might expect in terms of the 
relative duration of the blue phase, hence of the B/R ratio, we 
compare in the panels of fig.[2]the main physical and chemi- 
cal properties of the three models. In the left-bottom panel we 
report the variation with time of the effective temperature: we 
see that the drop of T e ff, indicating the return of the track to- 
wards the red, is postponed in the NACRE model, so that the 
blue phase lasts ~ 8Myr, to be compared to the ~ 6.5Myr dura- 
tion of the corresponding phase in the CF88 model, and to the 
~ 7.5Myr of the Kunz model. The ratio between the times spent 
in the blue and the red region of the diagram is consequently 
increase to tbi ue /t re d ~ 90% ( ~ 85% in the Kunz model). 

The longer duration of the blue phase is also confirmed 
by the three panels on the right, which show that the drop of 
CNO luminosity which is due to the general expansion of the 
structure, the phase of the rapid increase of the central density 
following helium exhaustion, and the increase of the central 
temperature are all slightly delayed in the NACRE model. The 
upper and middle panels on the left show the evolution of the 
chemical abundances of helium and carbon in the centre of the 
star. We note that the rate with which helium is consumed is 
higher in the CF88 model, while the carbon residual, as ex- 
pected is larger. The Kunz models has an intermediate behav- 
ior. 

To understand the reason of these differences, we focused 
our attention on the physical conditions of the central regions, 
particularly for what concerns the thermodynamical properties 
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Fig. 16. Upper: The variation with central temperature of the 
coefficients for nuclear energy release during helium burning 
within the interiors of two models presented in fig. ^] The 
lower lines refer to the energy due to the 12 C + a — > 16 O re- 
action, the middle lines to the reactions 3a and the upper lines 
represent the sum of the contributions, giving the global en- 
ergy release. Lower: Variation with central temperature of the 
central abundances of 4 He and I2 C within the same models. 

and the contribution of the two main reactions to the global 
nuclear energy release. 

First, we report in the two panels of fig. [21 the variation 
of central temperature and density in the very early phases of 
the He-burning phase (note the scale on the abscissa) as a func- 
tion of the central abundance of helium. For clarity reasons, we 
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limited our analysis to the NACRE and CF88 models, the Kunz 
model showing a behavior very similar to the NACRE. In both 
models we note a rapid increase of temperature followed by a 
temporary drop, which corresponds to a re-adjustment of the 
central regions following the rapid formation of a convective 
core of approximately M core ~ 0.2M G . The larger values of the 
NACRE crt2 C+a make the same demand for energy to be sat- 
isfied more easily, so that the temperature for a given helium 
abundance is slightly larger in the CF88 model: the differences 
in temperature in these early phases are very small, because 
the contribution of the quoted reaction to the global energy re- 
lease is very poor, but they will increase as helium burning pro- 
ceeds. The lower panel shows that no meaningfull differences 
are found in terms of density. On the basis of this, we expect a 
larger energy flux from 3a reactions in the CF88 model for a 
fixed helium abundance. 

In the upper panel of fig. [E|we show the variation with 
central temperature T c of the coefficient for nuclear energy 
generation e nuc i, together with the separate contributions from 
the 3a and the 12 C + a reactions; the lower panels shows the 
evolution of the chemical abundances of helium and carbon. 
To identify correctly the various phases, we marked with ver- 
tical lines the beginning of the excursion to the blue (which 
starts in both models when log(Tc) = 8.15), the bluest point 
(log(Tc) — 8.16), and the return to the red, pointing the end of 
the blue phase (log(Tc) = 8.2). We note that the values of Tc 
corresponding to the main evolutive stages are approximately 
the same for both models; therefore, it is the temporal varia- 
tion of Tc which determines the relative duration of the various 
phases, and the differences between the two models which we 
examine here. 

We see (fig. H6t that in both models &, a reaches a maximum 
when log{Tc) ~ 8.16, after which the helium depletion dom- 
inates on the temperature increase (we recall that e$ a ~ F 3 ). 
We note a larger 3a rate in the CF88 model which can be 
understood according to the results shown in the upper panel 
of fig. the CF88 model needs larger temperatures for a 
given helium abundance or, equivalently, as we may also see 
in the bottom panel of fig. ^] the same temperature corre- 
sponds to larger helium abundances, thus to larger 3a rates. 
The abundance of carbon, as expected, is seen in the lower 
panel to increase faster in the CF88 model, since less carbon 
is burnt via the 12 C + a — > 16 O reaction: yet, the energy re- 
lease due to this latter reaction is larger (see the upper panel) 
in the NACRE model, because of the larger cross-section pro- 
vided by the NACRE release. In both models the energy con- 
tribution due to 12 C + a burning is seen to increase for the 
whole duration of the He-burning phase, the temperature in- 
crease more than counterbalancing helium depletion (the de- 
pendence of er-c+a on X( 4 He) is weaker with respect to 3a re- 
actions, being eac+ a °c X( 4 He)X( l2 C)); when log(Tc) ~ 8.22, 
the energy release due to 12 C + a reaction exceeds the contri- 
bution provided by 3 a reactions. 

The comparison between the global energy release (see the 
upper panel of fig. 116b shows that for the same central tem- 
perature the CF88 model attains more nuclear energy by the 
two major reactions involving a particles; this determines a 
more rapid increase of the central temperatures, speeding-up 



the whole evolution. This means that the same central temper- 
ature is attained in the NACRE model somewhat later, so that 
the permanence in the blue region, which is terminated for both 
models when log(T c ) ~ 8.2, is slightly longer. 

8. Conclusions 

In this paper we present a detailed study of the helium burning 
phase within intermediate mass stars, and we apply our results 
to the LMC cluster NGC1866: we selected this cluster because 
it shows up a well populated clump region, allowing a statis- 
tically significant determination of the colour distribution of 
stars in the clump. The observational scenario indicates that 
the blue region of the clump is as much populated as the red, 
with an observed B/R ratio which is approximately 1. 

We tested the influence on the results obtained in terms of 
expected B/R ratio of some macro and micro-physics inputs 
adopted to calculate the models, namely: the scheme adopted 
for mixing within nuclear burning regions, the e-folding decay 
distance of convective velocities out of the formal border fixed 
by the Schwarzschild criterion, the description of convection, 
and the rate of the reaction 12 C + a — » 16 O. 

Our main findings are the followings: 

1 . The way nuclear burning is coupled to convective mixing 
within the central regions during the major phases of nu- 
clear burning has a strong influence on the results obtained, 
particularly in terms of the relative duration of the various 
evolutive phases. The models calculated with an instanta- 
neous mixing approximation burn helium very rapidly, so 
that the permanence in the blue region following the over- 
all contraction is shortened by the subsequent contraction 
of the nucleus, due to the lack of fuel; on the contrary, in the 
diffusive models helium burning is slower, so that a longer 
staying in the blue region of the clump is allowed, in better 
agreement with the observational scenario. 

The differences in terms of the relative duration of the blue 
phase amount approximately to a factor of ~ 2. Thus, our 
study of NGC1866 suggests that, during core He burning, 
mixing in a region of overshooting beyond the convective 
core boundary established by the Schwarzschild criterion 
should not be treated as instantaneous, but with some kind 
of decaying efficiency. 

2. The influence of ( , the parameter determining the e-folding 
distance of velocity decay within the radiatively stable re- 
gions, has been carefully tested. Larger values of £ al- 
low longer blue phases; models calculated with a diffusive 
scheme for mixing and £ = 0.03 lead to an expected B/R 
ratio of the order of ~ 75%. 

3. The cross-section of the reaction I2 C + a — > 16 O is the sec- 
ond major source of energy release during helium burning, 
so it may change the results obtained in terms of the dura- 
tion of various evolutive phases during helium burning. We 
find that turning to the NACRE (Angulo et al. 1999) value 
(which is approximately double than the value given by the 
Caughlan & Fowler (1988) prescription) leads to increase 
the permanence in the blue region, so that the expected per- 
centage of blue stars with respect to the red ones exceeds 
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90%, which is in excellent agreement with the observa- 
tional scenario. An analogous result is obtained by adopting 
the cross-sections for the quoted reaction given by Kunz et 
al. (2002) . 

4. The results obtained with the FST and MLT models for tur- 
bulent convection differ only in terms of the colour and the 
slope of the giant branch, but not for that concerning the ex- 
pected B/R ratio. No clear indication regarding convection 
can be got from this study. 
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